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1 .  Introduction 


In  general,  an  optical  limiter  keeps  the  power,  intensity,  energy  or  energy  density 
transmitted  by  an  optical  system  below  a  predetermined  maximum  value  that  is 
independent  of  the  size  of  the  input  pulse  while  maintaining  a  high  transmittance  at  low 
input  power.  The  many  applications  of  the  device  include  laser  power  regulation,  laser 
mode-locking,  optical  pulse  shaping,  signal  level  processing,  and  sensor/detector 
protection.  Important  applications  to  the  Air  Force  include  protection  of  sensors,  such  as 
human  eyes  of  pilots  and  photodetectors,  from  laser  damage.  In  order  to  determine  the 
parameters  required  for  optimization  of  optical  limiters,  we  are  developing  an  extremely 
sophisticated  numerical  code. 

Because  of  its  modularity  and  sophisticated  algorithms  this  code  has  wide  ranging 
future  applications  in  medical  optics,  nanooptics,  signal  processing,  multimedia, 
infomatics,  medical  imaging,  laser  damage,  photonic  crystals  and  industrial  devices. 
Therefore  it  has  potential  uses  in  both  the  military  and  civilian  domains.  Additionally, 
the  unique  ways  in  which  the  code  is  modularized  and  the  algorithms  performed  have 
important  Intellectual  Property  consequences.  We  believe  that  the  methods  used  here  are 
unique.  The  potential  applications  of  the  numerical  code  could  provide  numerous  new 
devices.  This  factor  is  significant  for  both  the  military  and  civilian  domains. 

Optical  limiting  now  requires  a  vast  set  of  properties  including  various  nonlinear 
materials,  different  physical  processes,  various  device  designs  including  a  number  of 
lenses  and  other  optical  devices,  optical  pulse  shaping,  and  the  characteristics  of  the  laser 
source  including  wavelength,  pulse  duration  and  pulse  repetition  rate.  This  range  of 
materials,  processes,  and  devices  involve  complex  optical  propagation  requiring  a  large 
diversity  of  parameters  that  cannot  be  fully  analyzed  without  the  use  of  various  numerical 
codes.  Along  with  new  materials  and  device  structures,  new  numerical  codes  will  be  key 
elements  in  the  development  of  optical  limiters.  Numerical  investigation  will  be  required 
to  make  additional  advances  by  examining  the  various  combinations  of  material, 
physical,  and  device  parameters  that  would  be  too  time  consuming  and/or  costly  to 
perform  in  the  laboratory. 

Therefore  one  of  the  key  features  to  the  successful  design  of  future  optical  limiting 
systems  will  be  sophisticated  numerical  codes.  As  stated  previously,  the  increasing 
complexity  of  the  materials,  device  geometry,  and  optical  beam  shaping  and  signal 
processing  necessitates  the  use  of  numerical  simulation  for  system  optimization.  The 
number  of  variable  parameters  that  will  have  to  be  optimized  can  approach  fifty  (50) 
different  variables.  Clearly  it  is  not  time  nor  cost  efficient  to  perform  all  of  these 
variations  experimentally.  The  numerical  codes  will  also  further  guide  quantum 
mechanical  calculations  and  molecular  engineering  in  order  to  optimize  the  possible 
molecular  structure  of  the  nonlinear  material.  This  code  is  modular  so  that  the  functions 
can  be  reused  or  eliminated  as  desired.  This  can  be  a  significant  feature  in  such  a  highly 
complex  code.  Furthermore  it  should  be  possible  to  use  additional  visualization  tools 
such  as  animation  that  could  enhance  the  understanding  of  the  optical  limiting  process 
and  enable  novel  advances  in  the  science  and  technology. 

A  second  area  of  investigation  involves  new  materials.  Accelerated  growth  in 
nanotechnology  has  led  to  materials  with  unique  physical  and  optical  properties  that  have 
never  been  observed  before.  These  developments  have  given  rise  to  potential 
applications  and  devices  with  physical  properties  that  can  far  exceed  those  of  traditional 
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materials/systems.  Advances  in  colloidal  chemistry  have  enabled  a  technology  to  create 
large  numbers  of  semiconductor  quantum  dots.  The  semiconductor  dots  can  be  deposited 
in  close-packed  films  using  organometallic  synthesis.  Recently,  large-scale  three- 
dimensional  arrays  of  semiconductor  dots  in  an  organic  medium  have  been  synthesized 
and  used  in  electron  transport  studies  of  films.  These  advancements  in  synthesis 
techniques  hold  promise  for  creating  other  classes  of  devices  with  different  electrical  and 
nonlinear  optical  properties. 

2.  Description  of  the  Methodology 

Because  of  the  complexity  of  the  beam  propagation  in  optical  limiting  systems 
numerical  methods  must  be  used.  The  numerical  methods  for  solving  the  optical 
propagation  along  the  path  from  the  incident  beam  to  the  detector  can  be  divided  into  the 
parts  shown  in  Fig.  1 . 


Figure  1.  Schematic  diagram  of  the  system. 


This  figure  demonstrates  the  simplest  optical  limiting  device  that  uses  a  single  nonlinear 
element;  however,  in  many  real  situations  additional  lenses,  limiting  materials,  and 
apertures  (not  include  in  this  figure)  could  be  included.  Region  I  shows  the  incident 
beam.  Regions  II  and  VI  are  thin  lenses  where  the  field  propagation  can  be  calculated 
using  a  simple  phase  shift.  After  propagating  through  the  first  lens,  the  beam  is  still 
Gaussian  shaped.  Therefore  in  Region  III  the  only  difference  in  the  field  is  the  beam 
waist.  Region  IV  contains  the  nonlinear  material  through  which  the  propagation  of  the 
field  may  include  nonlinear  absorption,  nonlinear  refraction,  and  dispersion.  However, 
after  propagation  through  the  nonlinear  medium  the  emerging  beam  is  usually  no  longer 
Gaussian  shaped  and  the  Huygen-Fresnel  propagation  formalism  is  required  for  beam 
propagation  in  regions  V  and  VII. 


2.1.  Regions  I-III 

The  paraxial  wave  equation  of  Gaussian  beam  can  be  written  as 
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A(r,r,z')  = 


^ Oe _ g'fr2  l2R(:!')e-r2/w2{z')e-rt  It l 

yll  +  Z,2/z20 


<j){z')  =  tan  1  (z'/z0) 


R(z')  =  z'  + 


where  z0  = 


A 


w(z')  =  w^\  +  z'2 1  zl 

is  Rayleigh  range,  also  known  as  confocal  parameter,  z'  is  defined  as 


the  distance  to  the  Gaussian  beam  waist  in  the  beam  propagation  direction,  w0  is  the 
radial  beam  waist ,  r  is  the  transverse  coordinate. 

A  Gaussian  beam  passing  through  a  lens  maintains  its  Gaussian  shape  because  the 
intensity  distribution  is  unchanged  by  a  thin  lens.  By  considering  the  incident  Gaussian 
beam  as  a  spherical  wave,  the  output  can  be  obtained  by  solving  the  ray  matrix  using  the 
ABCD  law.  After  the  lens,  the  emergent  Gaussian  beam  will  have  a  new  radial  beam 
waist  and  location  given  by 


.  / 

1  +  f2/zl 

,_*/  1 

™0>/l  +  f2/z20  ' 


where  d  is  the  distance  to  the  beam  waist  after  the  lens.  This  equation  describes 
focusing  a  collimated  beam.  Written  in  terms  of  the  f  (f  =  f/2wo),  Eq  (2)  becomes 

rf—  2^# 

1  +  (2h'0/#)2  /  Zq 

(3) 


2  Af#  1 

*  ^+(2wjy/z2 


The  f-number  is  often  f/5-f710.  In  Fig.  1,  d  is  the  distance  from  the  first  lens  to  the  focal 
plane  in  the  medium.  Therefore  the  focal  plane  can  be  adjusted  from  the  incident  plane 
to  the  exit  plane  in  order  to  best  utilize  the  nonlinear  material.  The  beam  waist  w'0  is 

assumed  to  be  a  collimated  beam  that  could  be  incident  from  far  away.  The  choice  of  the 
placement  of  the  focal  plane  could  also  be  significant  in  multiple  nonlinear  element 
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designs.  In  some  cases,  it  may  be  beneficial  to  have  multiple  elements  to  provide  optical 
limiting  for  multiple  wavelengths  and/or  multiple  pulses. 

2.2  Region  IV 

This  is  a  single  or  multiple  element  nonlinear  material  that  performs  the  functions  of 
beam  shaping  and  nonlinear  absorption.  In  the  nonlinear  material,  the  propagation 
operator  along  the  z-axis  is  written  as 


A(r,  t,z  +  A z)  =  exp  ^-iAzZ  ( r ,  t,  z)  J  A(r,  r,  z) 
d2  1  d 


(4) 


Z(r  ,t,z)  = 


2k 


—  ikaL  +  ko  Xnl  0,  r,z)  +  k  k2 


.  - L  '  '"0/tMA'  9  -9-/  ■  ••'•*'2  ^  2 

or  r  or  or 


A  unique  feature  of  our  code  is  that  it  discretizes  the  propagation  in  both  the  temporal  and 
spatial  domains  so  that  all  temporal  and  spatial  effects  can  be  observed.  In  this  respect  it 
is  the  most  detailed  code  to  date.  Additionally,  the  RSA  material  was  treated  as  a  five- 
level  model,  instead  of  the  three-level  model  approximation  used  in  previous 
calculations.  The  linear  operator  is  given  by 


The  nonlinear  operator  is  given  by 


2k 


i_i 

dr2  r  dr 
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The  propagation  equations  are  discretized  in  radial,  temporal  and  propagation  directions. 
This  method  enables  determination  of  all  spatial  and  temporal  features  of  the  beam 
propagation  through  the  entire  optical  system.  Thus  additional  features  may  be  observed 
with  this  code  that  are  not  possible  with  other  codes. 

2.3.  Regions  V  and  VII 

According  to  the  Huygens-Frensel  Principle  (HF),  every  point  of  the  wavefront  can  be 
considered  a  new  wave  source  that  mutually  interferes  with  one  another.  Using  the 
paraxial  approximation  and  the  HF,  the  field  can  be  expressed  as  a  plane  wave  given  by 


A  (z  +  Az,  r ,  t)  =  -7^—  exp 


ikr2 


\2lsz  j  o 


jj 

^r'dr  'A(z,r  ',r)exp 


ikr ' 
2A z 


krr' 

Az 


(7) 


where  Jn  is  zero  order  Bessel  function  and  Az  does  not  need  to  be  small. 


2.4  Region  VI 
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For  the  thin  lens,  in  the  paraxial  approximation  the  only  change  to  the  field  is  the  phase. 
The  field  after  the  lens  is  simply  expressed  as 


f 

A '( z ,  r,  r)  =  A(z,  r,  r)  exp 

V 


ikr2  N 

Jf; 


(8) 


where  /  is  the  focal  length. 


3.  Effects  of  Optical  Pulse  Trains 

To  date  numerical  codes  consider  only  a  single  incident  pulse.  Yet  most  lasers  emit 
trains  of  pulses  and  cw  lasers  can  be  externally  modulated  to  create  optical  pulse  trains. 
Therefore  it  is  of  interest  to  consider  the  effect  of  pulse  trains  on  optical  limiters.  A 
schematic  diagram  of  a  train  of  pulses  incident  to  an  optical  limiter  is  shown  in  Fig.  2. 


jin 


(a) 


^jutn 


Modulator 


(b) 


Figure  2.  Schematic  diagram  of  pulse  trains  incident  on  an  optical  limiter,  (a)  regularly 
modulated  pulse  train,  (b)  randomly  modulated  pulse  train. 


The  equation  for  the  incident  pulse  train  is  given  by 

N  r-(r-£r*)2 

"0VAF  - 2 - 

k= 1 


A(0,T)  =  ^A,cxp 


(9) 


where  rR  is  the  repetition  rate.  If  an  external  modulator  is  used,  then  the  pulse  train  can 
be  randomly  modulated  as  shown  in  Fig.  22  (b).  The  equation  for  this  case  is  given  by 


N 

A(0,T)  =  '£JAobkexp 


k=\ 


-{r-kTzf 


(10) 


where  bk  =0  or  1  depending  on  whether  the  kth  pulse  slot  contains  a  pulse  or  not. 

3.1  Investigation  of  Multiple  Pulses 
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The  use  and  significance  of  multiple  pulses  is  increasing  in  many  situations. 
In  addition  to  traditional  methods  such  as  mode-locking  or  Q-switching,  other  techniques 
are  possible  for  creating  pulse  trains  with  a  wide  range  of  repetition  rates  and  pulse 
widths  using  optical  modulators,  beam  splitters,  optical  delay  lines  and  booster 
amplifiers.  Our  model  and  numerical  method  is  applicable  to  a  wide  range  of  excitation 
conditions  however,  here  we  have  investigated  pulse  trains  from  5  to  20  pulses  at  two 
different  pulse  widths;  namely,  40  ps  and  1  ns  (1/e  half  width).  The  two  different  pulse 
widths  were  selected  in  order  to  demonstrate  and  validate  the  pulse  width  dependence  on 
the  intersystem  crossing  rate  of  the  molecular  system.  Since  we  found  that  there  was 
little  difference  in  the  general  physical  results  when  we  varied  the  number  of  pulses,  we 
present  results  for  five  pulses  as  examples. 

The  physical  values  used  for  our  calculations  are  taken  from  the  literature  for  the 
class  of  phthalocyanine  molecules  considered:  aoi=2.4xlO'18  cm2,  c?i2=3.0xl0'17  cm2, 
034=4. 8xl0'17  cm2,  kio=0.144  ns*1,  k2i=1.0  ps'1,  ki3=77.8  ps'1,  k3o=50.0  ms'1,  k43=1.0  ps'1 
The  other  parameters  used  in  our  calculation  are  Nt  (0.0048  nm'3),  beam  waist  (40  //  m), 
wavelength  (532  nm)  and  material  thickness  (0.6  mm). 

The  damage  threshold  for  various  photo  detectors,  including  the  cornea,  can  exhibit 
repetition  rate  dependence.  To  protect  such  detectors,  optical  limiters  will  need  to  be 
developed  that  are  operable  at  high  repetition  rates  or  when  subjected  to  complex 
excitations.  Recent  Z-scan  experimental  investigations  have  shown  that  multiple  pulse 
excitation  of  some  of  the  materials  being  investigated  for  optical  limiters  display 
dramatically  different  behavior  than  single  pulse  excitation.  As  lasers  become  more 
intense  and  exhibit  increased  repetition  rates,  a  detailed  quantitative  description  is 
required  to  describe  the  transient  response  of  the  material.  Fig.  3  shows  the  input  intensity 
normalized  to  Iq  as  a  function  of  time  normalized  to  t0. 


Normalized  Time  Domain  Plot 


Fig.  3.  Normalized  pulse  intensity  as  a  function  of  time  at  the  input  to  the  material. 
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Fig.  4.  Transmission  curve  as  a  function  of  input  energy  into  the  phthalocyanine 
chromophore  for  five  pulses  with  pulse  half-widths  of  1  ns  (solid  line)  and  40  ps  (dashed 
line),  pulse  separation,  tr  =5  and  L=0.6  mm. 


Fig  4  shows  the  transmission  curve  vs.  input  energy  for  both  the  1  nsec  pulse  sequence 
and  the  40  psec  pulse  sequence.  Transmission  is  defined  as  the  output  intensity 
integrated  over  r  and  t  divided  by  the  input  intensity  integrated  over  r  and  t.  The 
transmission  decreases  with  increasing  energy  for  both  pulse  width  sequences  within  the 
range  of  0.01  /u  J  to  100  juJ  due  to  excited  state  absorption  (reverse  saturable  absorption  - 
RSA).  Beyond  100  ju  J,  the  40  psec  pulse  sequence  experiences  saturation  whereas  the 
nsec  pulse  sequence  does  not.  The  onset  of  saturation  around  100  jj  J  for  the  psec  case  is 
due  to  the  depletion  of  ground  state  population,  which  reduces  the  excited  state 
absorption.  Intersystem  crossing  experienced  by  the  nanosecond  pulses  enables  the 
continued  RSA.  Since  the  relaxation  time  from  level  (3)  to  level  (1)  is  about  0.02  ms, 
the  successive  transfer  of  population  from  level  (3)  to  level  (4)  is  possible  with  the  optical 
pumping  for  the  nanosecond  pulses.  Since  the  excited  state  absorption  cross-section  <r34  is 

larger  than  that  of  the  ground  state  absorption  cross-section  cr01 ,  more  light  is  absorbed  if 

the  population  has  accumulated  at  level  (3).  These  effects  contribute  to  the  nonlinear 
absorption  up  to  lmJ  for  the  nanosecond  pulses. 

3.2  Investigation  of  temporal  and  radial  domains 

Next  we  investigate  the  temporal  pulse  shapes.  Fig.  5  (a)  shows  the  output 
pulse  intensities  for  40  ps  and  1  ns  pulse  sequences  as  a  function  of  time  and  radius  for 
various  input  energies.  The  first  energy  shown  is  0.01  //  J,  which  we  can  see  from  Fig.  4 
corresponds  to  the  linear  absorption  region.  In  this  case  the  normalized  peak  intensity  for 
each  of  the  pulses  is  0.45  in  agreement  with  Fig.  4.  The  next  input  energy  is  3.16  //  J  and 
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both  the  40  ps  and  the  1  ns  pulse  sequences  show  similar  behavior,  which  is  the  onset  of 
RSA.  At  this  onset  of  RSA,  it  is  interesting  to  note  that  the  trailing  pulses  are  more 
absorbed  than  the  leading  pulses  because  of  the  accumulation  of  population  in  state  1  due 
to  the  leading  pulses.  This  cumulative  effect  results  from  the  long  decay  rate  of  level  1. 
Since  the  excited  state  absorption  cross  section  is  larger  than  that  of  the  ground  state, 
there  is  more  than  enough  population  in  the  excited  state  to  allow  the  trailing  pulses  to  be 
excited  to  the  next  level,  which  causes  more  and  more  absorption  of  energy  from  the 
trailing  pulses.  At  an  input  energy  of  100  pi  J  we  see  the  onset  of  saturation  for  the  40  ps 
pulse  sequence  and  at  1000  pi  J  the  transmission  returns  to  the  linear  absorption,  which  is 
in  agreement  with  Fig.  4.  However,  the  behavior  of  the  1  ns  pulse  sequence  is  different. 
It  does  not  saturate  and  we  continue  to  see  the  decay  of  the  trailing  pulses  even  at  1000 
pi  J  input. 

We  now  examine  the  behavior  of  the  intensity  as  a  function  of  radius  for  the 
same  input  energies  as  above.  This  is  shown  in  Fig.  5  (b)  where  we  plot  the  intensity  of 
the  central  pulse  after  it  exits  the  material  as  a  function  radius  for  various  input  energies. 
At  an  input  energy  of  0.01  pi  J  the  transmission  value  is  0.45  in  agreement  with  Fig.  4. 
As  the  input  energy  increases  (e.g.  32  pi  J)  the  intensity  decreases  and  the  pulse  profile 
begins  to  flatten  in  the  radial  domain.  However  at  100  pi  J,  the  40  ps  pulse  shows  a  small 
increase  at  its  center.  This  energy  value  corresponds  to  the  minimum  transmission  value 
for  the  40  ps  curve  in  Fig.  4.  Comparison  of  the  corresponding  curve  for  the  1  ns  shows 
that  the  intensity  does  not  exhibit  this  increase  in  the  center  of  the  pulse.  Then  at  1000 
pi  I  we  see  that  the  40  ps  pulse  transmission  corresponds  to  the  linear  absorption  value; 
while  the  1  ns  pulse  remains  low.  The  radial  domain  plots  clearly  demonstrate  the  effects 
of  saturation  shown  in  Fig.4.  Therefore  both  the  temporal  and  spatial  domains  contain 
significant  information  about  the  laser-matter  interaction. 
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Fig.  5  (a)  Normalized  pulse  intensity  as  a  function  of  time  at  the  output  of  the  material 
for  five  pulses  with  pulse  half-widths  of  1  ns  and  pulse  separation,  tr  =5  and  r=0.  (b). 
Normalized  pulse  intensity  as  a  function  of  radius  at  the  output  of  the  material  for  five 
pulses  with  pulse  half-widths  of  40  ps  and  1  ns  and  pulse  separation,  tr  =5  and  t=0. 


The  population  densities  vs.  time  and  radius  are  shown  in  Fig.  6  and  Fig.  7  for  the  40 
psec  excitation  sequence.  The  populations  of  levels  3  and  4  remain  negligible  and  are  not 
shown.  The  input  energy  corresponds  to  the  onset  of  saturation,  100  ju  J.  Fig.  5  shows 
that  the  population  density  of  state  1  accumulates  with  successive  pulses.  The  shape  of 
the  population  density  vs.  radius  for  state  2  becomes  distorted  with  successive  pulses  and 
this  is  shown  in  Fig.  7  (a).  Fig.  7  (b)  shows  a  detail  of  the  radial  distortion  for  the  center 
pulse  in  the  sequence. 
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Fig.  6  Population  densities  for  levels  No  and  Ni  as  functions  of  time  and  radius. 
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Fig.  7.  Population  densities  for  N2.  (a)  three  dimensional  plot  of  intensity  as  function  of 
radius  and  time,  (b)  plot  of  intensity  as  function  of  radius  for  the  center  pulse  in  the  pulse 
sequence. 

Our  analysis  demonstrates  that  changes  in  the  radial  domain  need  to  be  considered  in 
order  to  understand  and  characterize  the  laser  matter  interaction.  This  particular  aspect  of 
the  analysis  is  traditionally  omitted  from  other  numerical  procedures  yet;  we  show  that 
the  radial  changes  can  be  significant. 
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4.  Semiconductor  quantum  dot  lattices 

4.1  Description 

Recent  results  show  that  semiconductor  dot  three-dimensional  arrays  in  organic  media 
lead  to  enhance  nonlinear  third  order  optical  susceptibilities.  When  semiconductor 
quantum  dots  are  placed  in  a  three-dimensional  lattice,  the  excitons  in  nearest-neighbor 
quantum  dots  interact  with  each  other  due  to  multipole  interaction.  This  interaction 
allows  the  presence  of  a  "transfer  Wannier-Mott  exciton  "  propagating  through  the  lattice. 
The  presence  of  this  transfer  exciton  helps  to  increase  the  optical  nonlinearity  of  the 
single  quantum  dots  to  overcome  the  limitation  of  the  single  dot  size  effect.  When  such  a 
lattice  is  put  in  an  organic  host,  this  transfer  Wannier  exciton  is  coupled  with  the  Frenkel 
exciton  of  the  organic  medium  and  at  resonance  of  the  two  excitons,  a  hybrid  exciton 
appears  in  the  system.  The  mixed  state  has  a  large  exciton  radius,  large  oscillator  strength 
and  very  large  optical  nonlinearity.  The  lattice  configuration  allows  tuning  the  optical 
properties  and  resonance  regions  as  well  as  the  value  of  the  nonlinearity  by  controlling 
the  dot  radius,  the  spacing  between  the  dots,  and  the  material,  especially  the 
semiconductor  dot  band  gap,  used  for  the  hybrid  arrays. 

The  role  of  systems  with  large  optical  nonlinearity  is  growing  in  importance. 
Important  recent  developments  in  technology  have  created  faster  and  more  intense  lasers, 
from  the  visible  to  the  far  infrared  thus  enabling  multiphoton  processes  across  a  broad 
range  of  wavelengths  from  the  visible  (~0.4-0.6  pm)  to  the  short  (0.8-2  pm),  mid  (3-8 
pm)  and  far  (>12  pm)  infrared  spectral  regions.  Two  photon  processes  are  particularly 
useful  in  many  applications  including  fluorescence  imaging,  optical  data  storage,  micro 
fabrication,  photoconductors  and  photovoltaics,  markers  for  genomes  and  proteins, 
biological  and  medical  detectors,  optical  limiters,  biomimetric  electromagnetic  devices 
and  nanopatteming  of  inorganic/organic  materials.  Consequently,  we  investigated  the 
two-photon  absorption  properties  of  this  array  of  semiconductor  dots  in  organic  materials 
with  its  hybrid  exciton. 


4.2  Hybrid-Exciton 

The  property  of  this  hybrid  exciton  formed  by  combining  the  “localized”  Frenkel 
exciton  with  a  “delocalized”  Wannier-Mott  exciton  has  been  of  considerable  recent 
interest.  The  Frenkel  exciton  is  the  excited  electronic  state  with  the  electron  and  the  hole 
situated  in  the  same  molecule  or  atom.  The  Frenkel  excitons  have  small  radii  and  very 
large  oscillator  strength.  In  constrast,  being  a  pair  of  electron  and  hole  in  the 
semiconductors,  the  large-radius  Wannier-Mott  exciton  in  semiconductors  is  relatively 
weakly  bound  due  to  typically  large  electron-hole  distance.  The  interaction  between 
Wannier-Mott  excitons  is  very  important  while  their  oscillator  strength  is  rather  weak. 
The  Wannier-Mott  exciton  wavefimctions  overlap  each  other  to  get  the  exciton 
resonance,  so  the  optical  nonlinearity  of  Wannier-Mott  exciton  can  be  large  at  rather  low 
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density.  These  differences  between  the  Frenkel  excitons  in  organic  materials  and  the 
Wannier  excitons  in  semiconductors  suggest  forming  systems,  which  have  a 
hybridization  of  both  the  excitons.  By  combining  an  organic  material  and  a 
semiconductor  in  one  heterostructure  system,  an  excitation  state  is  expected  to  have 
major  properties  of  both  kind  of  excitons  and  can  overcome  the  limitation  of  each  kind  of 
excitons.  Quantum  dots  have  attracted  a  lot  of  attention  because  of  their  three- 
dimensional  confinement.  The  electronic  structure  and  optical  properties  of  quantum  dots 
depend  sensitively  on  the  size  of  these  zero-dimensional  nanocrystals.  As  a  result  of  the 
confinement  effect  on  the  exciton  energy  levels,  the  oscillator  strength  of  optical 
transitions  in  quantum  dots  increases  with  the  size  of  crystallites  and  is  very  large  in 
comparison  with  the  bulk,  leading  to  an  increase  of  the  optical  nonlinearity. 

For  an  ideal  lattice  of  same-sized  quantum  dots  placed  in  an  organic  medium,  the 
total  Hamiltonian  is  given  by: 

H  =YJE:  (* K>„,  +  £  Erm  ( k)bl  btm  + 

n!  km  (11) 

Y  gnlm  0)0«>fe«  +  anlbL)  +  Y  *««'//'(£)(  «>»'/'  +  h.c .) 

nlmk  nn  'll ' 


where  gnim(k)  is  the  coupling  constant  of  Wannier  exciton  and  Frenkel  excitons,  and 
tm'i r(k)  is  the  hopping  constant  between  excitons  in  nearest  dots,  a*, ,  anl  ( b+km ,  bkm )  are 
creation  and  annihilation  of  Wannier  (Frenkel)  excitons,  respectively,  l  labels  the  exciton 
states  and  n  =  {i,j,h}  label  the  sites  in  the  quantum  dot  lattice,  E%  and  EF  are  the 
energies  of  the  Wannier  excitons  in  the  dots  and  the  Frenkel  exciton  in  the  medium, 
respectively. 

The  Hamiltonian  (Eq.  (1))  consists  of  the  free  Wannier  and  Frenkel  exciton  term,  the 
Wannier-Frenkel  exciton  coupling  term  and  hopping  term  between  the  Wannier  excitons 
in  nearest  dots.  Due  to  confinement,  the  Wannier  exciton  energy  in  the  dot  has  discrete 
value  according  to  the  zeros  of  the  Bessel  functions. 

The  hybrid  exciton  has  the  energy: 

E{k)  =  \{EF(k)  +  El  (*)}  ±  \ {[Ef  (k)  -  El  (k)f  +  4 G2 (/ k)}v 2  (12) 

and  the  wavefunction: 

*(*) = ui  mF  ml  ( k ) + vv{k)fw  (*)  03) 

where  vF;r(&),'P/F(k)  are  excited  states,  /^(0),/F(0)  are  ground  states  of  the 
Wannier  and  Frenkel  excitons  and  ut(k)  and  vr (k)  are  Bogoliubov  transformation 
coefficients.  The  Wannier-Frenkel  exciton  coupling  constant 
G(k )  =<  F(k)  |  Hint  |  W(k)  > ,  where  Hiat  is  the  coupling  of  the  field  created  in  the 
organic  medium  by  the  exciton  in  the  quantum  dot  with  the  transition  polarization 
operator  of  the  organic  medium,  and  the  hopping  constant  t(k)  =<  Wt  (k)  \  H d_d  \  W.  ( k )  > 

where  Hd_d  is  the  dipole-dipole  interaction  between  two  dipole  moments  in  the  two 
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nearest  dots  and  Wt  ( k )  is  the  wavefunction  of  the  exciton  in  the  dot  have  been  calculated 
in  detail.  The  energy  and  the  optical  properties  of  the  hybrid  exciton  depend  strongly  on 
the  values  of  G(k )  and  t(k). 

In  the  region  of  mixing,  both  the  Wannier  and  Frenkel  excitons  contribute  to  the 
oscillator  strength  of  the  hybrid  exciton.  And  one  of  the  strong  points  of  this  model  is  its 
ablity  to  tune  the  resonance  by  changing  the  materials  and  the  parameters  of  the  system  to 
obtain  the  desired  resonance  region  value  and  non-linearity  values.  For  example,  by 
choosing  a  semiconductor  with  a  band  gap  in  the  region  of  2.5  eV,  large  optical  third 
order  non-linearities  can  be  obtained  at  wavelengths  in  the  visible  and  near  IR  spectral 
region  where  many  lasers  exist,  such  as  the  Nd:YAG  laser  (1064  nm)  which  can  be 
frequency  doubled  (520  nm)  and  tripled  (355nm).  In  these  spectral  regions,  there  are 
many  applications  for  materials  with  large  optical  non-linearities,  such  as  optical 
communications  (e.g.  optical  switches)  and  laser  hardening  (e.g.  optical  filters  and  optical 
limiters).  Both  the  real  and  imaginary  parts  of  the  non-linearity  are  potentially  of  interest. 
The  real  part  gives  rise  to  self-focusing  or  self-defocusing  depending  on  the  sign  of  the 
non-linearity;  while,  the  imaginary  part  gives  rise  to  two-photon  absorption.  Figs.  8  and 
9  show  schematic  diagrams  of  thin  film  structures  of  inorganic-organic  materials.  Fig.  8 
shows  one  type  of  semiconductor  in  an  organic  medium  and  Fig.  9shows  a  multi-layered 
structure  in  which  different  semiconductors  are  used.  For  example,  one  layer  may  be  a 
self-focusing  layer  and  the  other  layer  may  be  a  two-photon  absorbing  layer.  Various 
materials  and  structures  may  be  used  for  different  applications. 
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Fig.  8.  Schematic  diagram  of  a  thin  film 
organic-inorganic  structure. 
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Fig.  9.  Schematic  diagram  of  a  multi-layered  thin  film  structure  with  different  types  of 
semiconductors  and  different  dot  spacing  in  each  layer. _ 


Due  to  the  confinement  effect  of  the  exciton  as  well  as  the  transfer  exciton 
coupling  between  dots  in  the  array,  the  Wannier  transfer  exciton  may  achieve  the 
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oscillator  strength  comparable  to  that  of  Frenkel  exciton,  and  so  the  hybrid  exciton  has 
very  large  oscillator  strength.  In  the  presence  of  an  optical  field,  the  third  order 
susceptibility  has  been  calculated  using  the  standard  perturbation  theory.  The  increase  of 
the  oscillator  strength  as  well  as  the  increase  of  coherence  length  leads  to  a  large  figure 
merit  of  the  hybrid  exciton  and  the  large  nonlinearity.  The  imaginary  part  of  the  third 
order  nonlinearity  is  given  by 
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(14) 


where  lc  is  the  coherence  length,  y±  and  yD  are  the  transverse  and  longitudinal 

relaxation  constants  of  the  excitonic  transition,  respectively,  R  is  the  dot  radius,  d  is  the 
dot  to  dot  separation,  Vm  is  the  volume  of  the  organic  host,  Vc  is  the  volume  of  one  cell  in 
the  organic  lattice,  pp  is  the  optical  transition  dipole  moment  of  the  organic  molecule,  a  is 
the  Bohr  radius,  and  hnr  is  the  lowest  excitation  energy  of  the  hybrid  exciton.  The 
parameters  used  for  the  calculation  are  hex  =  25  meV  (i.e.  CdS),  a  =  30  A,  Vc  =  (5  A)3, 
Vm  =  (5000  A)3,  pF  =  5  D  (lD=3.335xlO'30  Coulomb.meter),  y±=5.0xl012  s*1  and 
yQ  ==2.27x1 09  s'1.  Our  calculations  assumed  CdS  quantum  dots  with  a  dot  size  of  about 
30-100  A  and  a  dot-to-dot  spacing  ranging  from  2  to  5  times  the  dot  size. 

5.  Application  of  Hybrid  Material 


As  one  example  of  an  application  of  these  hybrid  structures,  we  investigate  their 
potential  for  optical  power  limiting  or  optical  filtering.  Figure  10  shows  plots  of  the 
imaginary  part  of  the  third  order  nonlinear  susceptibility  as  a  function  of  wavelength. 
The  resonances  occur  near  500  nm  and  the  non-linearity  is  large  even  at  off-resonance 
wavelengths  including  the  near  IR  region.  Future  studies  will  investigate  the  real  part  of 
the  non-linearity  that  gives  rise  to  the  nonlinear  index  of  refraction  coefficient.  As  stated 
previously  self-focusing  or  self-defocusing  results  from  the  real  part  of  the  non-linearity. 
For  example,  very  strong  self-focusing  may  enable  smaller  spot  size  for  laser  induced 
micro  fabrication.  This  may  lead  to  improved  methods  for  making  nanostructures  such  as 
photonic  crystals. 
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Fig.  10.  Plots  of  Imaginary  part  of  the  third  order  nonlinear  susceptibility  as  a  function  of 
wavelength.  The  dotted  curves  correspond  to  R=  50  A  and  d=100A,  the  dash  dot  curves 
correspond  to  R=50  A  and  d=150  A,  and  the  solid  curves  correspond  to  R=40  A  and 

d=160  A. 
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There  are  four  resonances  that  occur  in  the  visible  spectral  region.  Fig.  10  shows  the 
imaginary  part  of  the  third  order  susceptibility  as  a  function  of  wavelength  for  three  ratios 
of  d/R;  namely,  d/R  =2,  3,  4.  Since  quantum  confinement  restricts  the  size  of  the  dot  to 
-3-10  nm,  values  between  4  nm  and  5  nm  were  selected  as  initial  examples.  The  curves 
show  that  as  the  ratio  d/R  increases  the  non-linearity  decreases.  At  off-resonance 
wavelengths,  our  calculations  show  that  about  an  order  of  magnitude  (i.e.  -  factor  of  10) 
increase  in  the  non-linearity  is  possible  when  the  d/R  ratio  is  reduced  by  only  a  factor  of 
2  (e.g.  from  4  to  2).  Thus  a  significant  change  in  the  non-linearity  can  be  achieved  by 
relatively  small  changes  in  this  ratio.  These  results  indicate  the  ability  to  significantly 
tune  the  non-linearity  by  changing  the  d/R  ratio.  Further  studies  are  required  to  examine 
the  effects  of  various  semiconductors  and  additional  changes  in  the  dot  size  and  dot-to- 
dot  distance. 

Devices  containing  optical  limiting  materials  are  often  used  for  laser  protection  in 
the  visible  and  infrared  spectral  regions  for  optical  (human  eyes)  and  optoelectronic 
detectors.  Many  different  types  of  materials  have  been  investigated  for  optical  limiters 
such  as  organics,  fullerenes  and  semiconductors.  With  the  recent  advancements  in  laser 
technology,  additional  materials  are  required  for  optical  limiters.  Typical  optical  limiters 
in  the  visible  and  near  infrared  regions  reduce  the  transmission  of  the  incident  laser  beam 
by  about  an  order  of  magnitude.  Yet,  for  many  applications,  it  is  advantageous  to  have  a 
single  element  that  can  reduce  the  transmitted  energy  by  several  orders  of  magnitude 
(-1000).  An  important  physical  mechanism  used  for  optical  limiting  is  two-photon 
absorption.  The  imaginary  part  of  the  susceptibility  gives  rise  to  the  coefficient  for  two- 
photon  absorption  (a2)  and  is  given  by 

(15) 

At  this  time,  our  calculation  of  the  two-photon  absorption  fa)  coefficient,  for  the 
inorganic-organic  material  described  in  this  paper,  shows  that  it  significantly  exceeds  (by 
two  to  five  orders  of  magnitude)  that  of  other  traditional  materials  in  the  visible/near 
infrared  spectral  regions.  Future  work  will  investigate  other  materials  for  comparison. 

For  example,  we  obtain  values  of  a2  -1 .2  xlO'9  m/W  {1}  for  CdS  quantum  dots  (R-100 
A,  d-200  A)  in  an  organic  medium  at  a  wavelength  of  800  nm.  This  value  exceeds  a2  for 
other  materials  in  the  visible  and  near  infrared  spectral  regions.  For  example,  a2  -10"14 
m/W  {II}  (532nm)  for  2, 2’-(9, 9-dihexyl)  bifluorene,  a2-6xl0"u  m/W  {III}  (800nm)  for 
green  fluorescence  protein,  and  a2~10'"m/W  {IV}  (1060nm)  for  trans-4-[p-(N-ethyl- 
N-hydroxethylamino)styrl]-N-methypyridinium  tetraphenylborate. 
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Table  I.  Values  of  the  two-photon  absorption  coefficient  as  a  function  of  wavelength 


Material 

a2  (m/W) 

Wavelength  (nm) 

{1} 

1.2x  10‘9 

800 

{11} 

X 

o 

532 

{III} 

6  x  10'11 

800 

{IV} 

lx  10'11 

1060 

We  calculated  the  transmission  through  a  hybrid  semiconductor  dot-organic  thin  film 
using  the  two-photon  absorption  parameters  obtained  from  Eq.  (15).  Potasek  et.  al.  have 
developed  theoretical  models  and  numerical  codes  for  laser  pulse  propagation  in 
nonlinear  absorbers.  Using  numerical  simulation  we  calculated  the  transmission  of  a 
Gaussian  laser  beam  through  a  hybrid  3D  semiconductor  array-organic  film  with  a 
thickness  of  about  500  pm  at  an  incident  wavelength  of  800  nm.  The  figure  shows  more 
than  three  orders  of  magnitude  (>1000)  reduction  in  transmission.  This  reduction  is 
greater  than  that  obtained  using  traditional  materials. 


Fig.  1 1  Transmission  as  function  of  input  energy 

6.  Conclusions 
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Considerable  new  efforts  were  added  to  the  investigation  of  nonlinear  absorbers  that 
are  critical  to  laser  protection.  New  advances  in  femtosecond  lasers  with  short  duration 
optical  pulses  and  fast  repetition  rates  increase  the  possibility  of  laser  damage  in  this  new 
regime.  Our  numerical  effort  has  expanded  into  these  new  areas. 

Additionally  we  show  that  an  array  of  semiconductor  quantum  dots  in  an  organic  host 
leads  to  large  two-photon  absorption.  The  optical  nonlinearity  depends  on  the 
semiconductor,  the  dot  size  and  the  dot-to-dot  spacing.  Using  numerical  simulations,  we 
demonstrate  that  large  optical  limiting  is  possible  using  thin  films  of  this  hybrid  material. 
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